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ABSTRACT 

We study, analytically and numerically, the energy input into dark matter mini-haloes 
by interactions with stars. We find that the fractional energy input in simulations 
of Plummer spheres agrees well with the impulse approximation for small and large 
impact parameters, with a rapid transition between these two regimes. Using the 
impulse approximation the fractional energy input at large impact parameters is fairly 
independent of the mass and density profile of the mini-halo, however low-mass mini- 
haloes experience a greater fractional energy input in close encounters. We formulate 
a fitting function which encodes these results and use it to estimate the disruption 
timescales of mini-haloes, taking into account the stellar velocity dispersion and mass 
distribution. For mini-haloes with mass M < O(10~^Mq) on typical orbits which pass 
through the disc, we find that the estimated disruption timescales are independent of 
mini-halo mass, and are of order the age of the Milky Way. For more massive mini- 
haloes the estimated disruption timescales increase rapidly with increasing mass. 

Key vifords: Dark matter. Galaxy-structure 



1 INTRODUCTION 



In Cold Dark Matter (CDM) cosmologies structure forms 
hierarchically; small haloes form first, with larger haloes 
forming via mergers and accretion. The internal structure 
of haloes is determined by the dynamical processes, for in- 
stance tidal stripping and dynamical friction, which act on 
the component sub-haloes. Numerical simulations find that 
substantial amounts of substructure survive within larger 
haloes (Klypin et al. 1999; Moore et al. 1999) with the num- 
ber density of sub-haloes increasing with decreasing mass, 
down to the resolution limit of the simulations. Natural ques- 
tions to ask are: what are the properties of the first dark 
matter (DM) haloes to form in the Universe, and do signif- 
icant numbers survive to the present day? 

Weakly Interacting Massive Particles (WIMPs) are one 
of the best motivated DM candidates. They generically have 
roughly the required present day density, and supersymme- 
try (an extension of the standard model of particle physics) 
provides a well-motivated WIMP candidate, the lightest 
neutralino (see e.g. Bertone, Hooper & Silk 2005). Numer- 
ous experiments are underway attempting to detect WIMPs 
either directly (via elastic scattering off target nuclei in the 
lab) or indirectly (via the products of their annihilation). In 
both cases the expected signals depend critically on the DM 
distribution on small scales; direct detection probes the DM 
on sub-milli-pc (AU) scales, while clumping may enhance 
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the indirect signals. Therefore the fate of the first dark mat- 
ter haloes, and the resulting dark matter distribution on 
small (sub-pc) scales, is important for practical reasons too. 

Studies of the microphysics of WIMPs show that kinetic 
decoupling and free-streaming combine to produce a cut-off 
in the density perturbation spectrum for generic WIMPs at 
a comoving wavenumber k ~ 0(1 pc~^) (Hofmann, Schwarz 
& Stocker 2001; Schwarz, Hofmann & Stocker 2001; Berezin- 
sky, Dokuchaev & Eroshenko 2003; Green, Hofmann & 
Schwarz 2004, 2005; Loeb & Zaldarriaga 2005; Berezinsky, 
Dokuchaev & Eroshenko 2006) which corresponds to a mass 
of order IO'^A/q. Analytic calculations, using linear theory 
and the spherical collapse model, find that the first typical 
(i.e. forming from 1-a fluctuations) haloes form at 2; ~ 60 
and have radius R ~ 0(0.01 pc) (Green et al. 2004, 2005). 

Diemand, Moore & Stadel (2005) carried out numer- 
ical simulations using as input the linear power spectrum 
for a generic WIMP with mass m^ = 100 GeV (Green et 
al. 2004, 2005). They used a multi-scale technique, twice re- 
simulating at higher resolution an 'average' region selected 
from a larger simulation. These simulations confirmed the 
analytic estimates of the mass and formation red-shift of 
the first mini-haloes and also provided further information 
about their properties, in particular the density profiles of 
sample mini-haloes. The simulations were stopped at z ~ 26 
when the high-resolution region began to merge with the 
lower resolution surroundings, and so the subsequent evolu- 
tion of the mini-haloes has to be studied separately. 

Extensive work has been done on the dynamical evolu- 
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tion of more massive (M > 10^ Mq) substructure (e.g. Zent- 
ner & Bullock 2003; Taylor & Babul 2004; Oguri & Lee 2004; 
Penarrubia & Benson 2005). The physics of mini-haloes is 
significantly different to these more massive haloes, however. 
Firstly, the first generation of mini-haloes form monolithi- 
cally, rather than hierarchically. Secondly, the amplitude of 
the density perturbations on these scales is a very weak func- 
tion of scale, so that haloes with a range of masses form at 
the same time. Finally, as well as being subject to the same 
dynamical processes as larger sub-haloes (e.g. tidal strip- 
ping, interactions between sub-haloes[j) mini-haloes can lose 
energy, and possibly be completely disrupted, via interac- 
tions with compact objects such as stars. Various authors 
have used the impulse approximation to investigate the dis- 
ruption of Af ^ 1O~^M0 mini-haloes due to encounters with 
stars (Diemand et al. 2005; Zhao et al. 2005a; Moore et al. 
2005; Zhao et al. 2005b; Berezinsky et al. 2006; Goerdt et 
al. 2006; Angus & Zhao 2006). The resuhs of these stud- 
ies range from most of the mini-haloes surviving disruption 
(Diemand et al. 2005; Moore et al. 2005) to most of the 
mini-haloes whose orbits pass through the solar neighbour- 
hood being destroyed (Zhao et al. 2005a). A definitive study 
will have to combine accurate calculations of the response 
of mini-halos to individual interactions with simulations of 
mini-halo orbits in a realistic Galactic potential, (see Moore 
(1993), Zhao et al. (2005b) for work in this direction). 

In this paper we use N-body simulations to investigate 
the accuracy of the impulse approximation for calculating 
the energy input into a mini-halo by an interaction with a 
star. We formulate a fitting function which matches the re- 
sults of the simulations and use it to estimate the timescales 
for one-off and multiple disruption as a function of mini-halo 
mass. We caution, and discuss in more detail below, that it 
is actually the mass loss which is key to determining the 
extent to which a mini-halo is disrupted. The relationship 
between the energy input and the mass lost in an interaction 
is a complex, and to some extent unresolved, problem (see 
e.g. Aguilar & White (1985) and Goerdt et al. (2006)). 



2 PREVIOUS CALCULATIONS 

The duration of a typical star-mini-halo encounter is far 
shorter than the dynamical time scale of the mini-halo, 
therefore the impulse approximation holds and the inter- 
action can be treated as instantaneous (Spitzer 1958). More 
specifically, the validity of the impulse approximation can 
be parameterised by the adiabatic parameter (Gnedin & 
Ostriker 1999) x = lot where ui = a{b)/b is the orbital fre- 
quency of particles at a distance b from the centre of the 
mini-halo and r = 2R/v is the duration of the encounter. 
The impulse approximation is valid if a; <g 1, or equivalently 
R/b <g v/a(b). As the typical relative velocity of encounters 
{v ~ 0(10 — 100 kms^^)) is far larger than the mini-halo 
velocity dispersion (cr(_R) ~ 0(1 ms~^)), then only for very 
rare, slow interactions at very small impact parameters will 
the impulse approximation be violated. The change in the 
velocity of a particle within an extended body of radius R 
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at position r relative to the centre of the body due to an 
impulsive interaction with a perturber of mass A/* , moving 
with relative velocity v at an impact parameter b, such that 
6 > i?, is given by (Spitzer 1958) 
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where ev and eb are unit vectors perpendicular to v and b 
respectively. The energy input, per unit mass, for an individ- 
ual particle is 5E = v.(5v) -|- 0.5((5v)'^, and the total energy 
input into the body is then found by integrating over the 
density distribution. For a spherically symmetric body the 
first term averages to zero and, using the approximation that 
(r.Gv)^ — (r.Gb)^ ~ f^/-^, the total energy input is given by 
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is the root mean square radius. 

For small impact parameter interactions, b/R -^ 0, (e.| 
Gerhard & Fall 1983, Carr & Sakellariadou 1999) 
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so that the energy input is given by 
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is the root mean square inverse radius. Carr & Sakellari- 
adou (1999) (drawing on Gerhard & Fall (1983)) interpolate 
between the b^ R and b <^ R regimes using 
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Moore (1993) simulated encounters between globu- 
lar clusters, modelled with King profiles, and massive (~ 
10* Af©) black holes. He found that the energy input was 
well fitted by 



AE{b) = 



AE{b = 0) 
[1 + (b/R)]' 



(9) 



This fitting function will, however, only reproduce the 
asymptotic limits, equations ((5} and ([SJ, if a^ « 9/3^/4, 
which need not be (and we will see is not) the case in gen- 
eral. A simple modification to equation ([SJ 
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halo 


profile 


a2R2 (p^2) 


p^/R^ (pc-2) 


E/M (crg/Me) 


M (Mq) 


'•p/c/s (pc) 


1 
1 
1 


plumincr 
CIS 

NFW 


2.6 X 10"'' 
3.2 X 10-4 
3.2 X lO-'' 


1.5 X 10* 
3.2 X 10* 


-1.5 X 10^* 
-7.2 X 10^'^ 
-7.8 X 10^^ 


1.0 X 10-* 
6.0 X lO--"' 
6.7 X lO--"' 


0.013 

0.0032 

0.0014 



2 


plumincr 


3.2 X lO--""' 


7.0 X 10* 


-2.1 X 10^7 


2.1 X lO-'' 


0.0092 


2 


CIS 


2.7 X lO--""' 


1.2 X 10^ 


-1.1 X IQS'^ 


1.6 X lO-'' 


0.0032 


2 


NFW 


2.7 X 10-"' 


- 


-5.1 X 10^'^ 


1.4 X 10-^ 


0.011 


3 


plummcr 


3.0 X 10-"' 


8.1 X 10* 


-1.0 X lO^^^ 


1.3 X 10-^ 


0.0076 


3 


CIS 


2.4 X 10-"' 


2.0 X 10^ 


-1.0 X lO^'i' 


1.2 X lO-'' 


0.0019 


3 


NFW 


2.6 X 10-"' 


- 


-5.2 X 1037 


1.4 X 10-6 


0.0076 



Table 1. Structure parameters, binding energy per unit mass and mass of the best-fit profiles for the 3 sample haloes (using only the data 
at radii greater than the force softening). For the energy and mass calculations a sharp cut-off is taken at ij = r20o{z = 26) = 0.03 pc 
and 0.008 pc for halo 1 and haloes 2 & 3 respectively. The final column is the Plummer, core or scale radius for the Plumincr, CIS and 
NFW profiles respectively. 
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[1 + (6A-i/4/i?)]* ' 

with A — 4a^/9/3^, produces a function with the correct 
asymptotic limits in general. 

For later convenience we write the fractional energy in- 
put, which is simply AE{b) divided by the total energy of 
the mini-halo E = 'yGM / R, where 7 is a constant of order 
one which depends on the density profile (and cut-off radius 
if one is imposed) , as 
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where {AE{b)/E)fid is the fractional energy input in an in- 
teraction with a fiducial star with mass M* = IMq and 
relative velocity v = 300 km s-^. 



3 APPLICATION TO MINI-HALOES 

3.1 Density profiles 

We use three benchmark density profiles: 



p{r 



(i) Plummer sphere 
Po 



[1 + ir/r^yf' 



(12) 



This profile (Plummer 1915), which has a central core and 
asymptotes to r~^ at large radii, is commonly used to model 
star clusters. It is not a good fit to simulated CDM halos 
or subhalos, however it is a convenient choice for testing the 
impulse approximation against numerical simulations as it 
has a simple form for the density and velocity distributions 
(see Aarseth, Henon & Weilan 1974). In addition, the rapid 
fall off of the density at large radii means that there are no 
subtleties involved in imposing a truncation radius (see e.g. 
Kazantzidis, Magorrian & Moore 2004), and the Plummer 
sphere is also stable when isolated, 
(ii) cored isothermal sphere (CIS) 
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The cored isothermal sphere is a better (although still not 
good) approximation to the mini-haloes, is amenable to an- 



alytic calculations and allows us to investigate the impact 
of a central core and more gradual fall-off at large radii to 
the fractional energy input. 

(iii) Navarro, Frenk, White (NFW) 



p(r) = 



Po 



(r/r,)[l-h(r/r,)]2 



(14) 



The Navarro, Frenk, White profile (Navarro, Frenk & White 
1996, 1997) fits the density distribution, outside the very 
central regions, of simulated galactic scale and larger dark 
matter haloes well and is often used to model massive dark 
matter haloes. However, mini-haloes form monolithically, 
rather than by hierarchical mergers like 'standard' dark mat- 
ter haloes, and it is not clear that they will have the same 
density profile. The NFW profile does, however, provide a 
reasonably good fit to the density profiles of the mini-haloes 
from Diemand et al.'s simulations. 

We find the best fit for each of these profiles for the 
three typical haloes in Fig. 2 of Diemand et al. (2005) using 
only the data points (density averaged within radial bins) 
at radii greater than the force resolution (using all the data 
points does not significantly change the best fit parameters). 
We refer to the haloes denoted by squares, stars and circles 
in their figure as haloes 1, 2 and 3 respectively. 

The CIS and NFW profiles have infinite mass and en- 
ergy unless a cut-off radius is imposed by hand. For definite- 
ness, and to allow comparison with previous work on mini- 
halo disruption, we use the radius at which the halo density 
is 200 times the cosmic mean density at z = 26 when Die- 
mand et al.'s simulations are stopped and the sample haloes 
studied. For halo 1 r2oo(z = 26) = 0.03 pc while for haloes 
2 and 3 r2oo{z = 26) — 0.008 pc. For the Plummer profile a 
cut-off is not in principle needed, however we use the same 
values for the radii for consistency. 

In Table 1 we give the values of the Plummer, core and 
scale radii (as appropriate) , the mass (M) , initial energy per 
unit mass (E/M) and the structure parameters, a^R^ and 
/3^/i?^, Our values of the structure parameters are slightly 
different to those of Carr & Sakellariadou (1999) as they de- 
fine the cluster radii differently. Haloes 2 and 3 have roughly 
the same mass ~ 10~^Mq depending at the tens of per-cent 
level on the profile used. Halo 1 is a factor of ~ 50 more 
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massive Q We calculate the total energy by calculating the 
potential, and hence the velocity dispersion and kinetic and 
potential energy densities, from the density profiles. The 
truncation at finite radii means that the resulting haloes 
are not in virial equilibrium. For all three profiles for halo 
1 and the CIS and NFW profiles for haloes 2 and 3 the de- 
viation is relatively small. The best fit Plummer spheres for 
haloes 2 and 3 have r^ ~ 0(^200(2 = 26)) and the resulting 
systems are far from virial equilibrium. 

The values of a^, which parameterises the energy in- 
put for large impact parameter encounters, only vary by a 
factor of '^ 2 between different haloes and density profiles 
refiecting the fact that the sample halos have similar mean 
densities. However /3^, which parameterises the energy in- 
put in the 6-^0 limit varies significantly and is in fact 
infinite for the NFW profile. It can be seen from the def- 
inition of /3^, equation ((6)1, that p'^ is formally infinite for 
any profile with a central cusp p{r) ~ r^^ with 7^1. The 
WIMP density can not in fact become arbitrarily high in the 
central regions of a mini-halo; if the density becomes suffi- 
ciently high the WIMPs will annihilate, reducing the density 
to some maximum value pmax so that the halo has a (small) 
core: p{r) = pmax for r < rcorc. The density and size of the 
core can be estimated by calculating the density for which 
the annihilation time scale is less than the Hubble time (c.f. 
Berezinsky, Gurevich & Zybin 1992) 



Pmax < ''"xx^ > 
2m^ 



IQiOyr 



(15) 



Using 'typical' values for the WIMP mass and velocity 
averaged cross-section, rrix ~ 100 GeV, < cfxx^ >= 3 x 
10"^^m^s"\ we find p^ax = 4 x 10"^^kgm"^ = 4 x 
10"pc(z = 0). For the best fit NFW profiles rcorc ~ 10"^'' pc. 
Taking this effect into account leads to finite values for /3^, 
but they are still large (~ 10 ). The energy input only 
reaches its asymptotic value, however, for tiny (~ rcorc), 
and hence extremely rare, impact parameters. 



3.2 Simulations of star — mini-halo encounters 

We use the DRAGON smooth particle hydrodynamics code 
(e.g. Goodwin, Whitworth & Ward- Thompson 2004a,b; 
Hubber, Goodwin & Whitworth 2006) with hydrodynamics 
turned-off as an A''-body code. DRAGON uses a Barnes-Hut 
(1986)-type tree and we set the opening angle to be small to 
increase the accuracy of the force calculations between DM 
particles. The forces between DM particles and the star are 
all computed by direct summation. This physical situation, 
interaction of an extended body with a far more massive 
compact object, has, to our knowledge, not been studied 
numerically before and we carried out extensive testing to 
ensure the reliability of the results. In particular, the masses 
of the DM particles are a factor ~ 0(10^) less massive than 
the perturbing star, requiring numerical care to be taken. 

We generate the initial conditions for the Plummer 
sphere mini-haloes using the prescription of Aarseth et al. 



(1974), assuming that the haloes are initially in virial equi- 
librium. Left isolated, the mini-haloes remain in equilibrium, 
and the energy conservation of the code is ~ 10~^ over 
timescales far in excess of a typical mini-halo-star interac- 
tion timescale (~ 50 kyr). A star of mass M* is then placed 
Ipc away from the halo approaching it at velocity v, with 
an impact parameter b. 

We conduct simulations with A*' = 5000 DM particles 
with a Plummer force softening between DM particles of 
e — 10~^ pc. The forces due to the star are softened with a 
significantly smaller softening length of 10~* pc. The soft- 
ening between DM particles is rather large, but we wish 
to subdue any 2-body interactions between DM particles. 
Tests conducted with e = 10~* and 10~^ pc show no differ- 
ence in the results. Similarly, increasing the particle num- 
bers to A'^ = 10000 and 20000 we found no significant (or 
systematic) changes. This convergence is not surprising as 
the energy input is entirely due to the encounter with the 
star whose force is accurately calculated with a low soften- 
ing length, and we are only concerned with the energy input 
to the halo, and not in the details of relaxation and/or mass 
loss after the impulse has occurred (which will involve inter- 
actions between the halo particles and may require a larger 
number particles for convergence e.g. Goerdt et al. 2006). 

We ran a large ensemble of simulations covering a wide 
range of Af* —v~b parameter space: 0.215 < M*/Mq < 30, 
1 < «/(lkm s-i) <40oE|and-5<logio(fe/lpc) < 1. With 
N = 5000, each simulation took an average of 20 minutes 
on a desktop PC. 



3.3 Fractional energy input 

The fractional energy input, {AE{b)/E)fid, in an interaction 
with a fiducial star with mass M* = IMq and relative veloc- 
ity V — 300kms~^ is plotted in Fig. 1 for the best fit plum- 
mer sphere for halo 1. This fiducial velocity was chosen as an 
isothermal sphere with circular velocity Vc = 220 km s~^ (i.e. 
representing the Milky Way) has root mean square speed of 
270 kms"'^. In reality, interactions will have a range of ve- 
locities and perturber masses. In Figs. [2] and |3] we plot the 
fractional energy input as a function of relative velocity and 
perturber mass, showing that it scales as v~'^ and M^ re- 
spectively as expected from equation ([2]). In Fig. [2] we also 
show that the v~^ proportionality is independent of the im- 
pact parameter and holds down to very small (0(1 km s~^)) 
relative velocities. 

For a given perturber mass and relative velocity, we see 
from Fig. 1 that the large and small b asymptotic limits are 
in excellent agreement with the full analytic calculation us- 
ing equation ((Sjl. As expected, the original fitting function 
significantly over-estimates the energy input at large b. The 
modified fitting function, designed to reproduce the asymp- 
totic limits, matches well the calculation using equation ([8} 
for all b. In the simulations, however, the transition between 
the b <^ R and b ^ R regimes happens very rapidly and the 
energy input is well approximated, for all 6, by the minimum 
of the asymptotic limits: 



^ It appears that the 5.1 x 10 ^Mq for halo one in the caption 
of fig. 2 of Dieniand et al. (2005) is a typo and should be 5.1 x 
lO-^M©. 



•^ Although interactions with relative speeds at the lower end of 
this range are extremely rare we consider them in order to test 
the validity of the impulse approximation. 
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Figure 1. The fractional energy input, {AE{b)/E)fi^, in an in- 
teraction with a fiducial star with mass Af* = IMq and relative 
speed V = SOOkms"^ for the best fit Plunimer profile for halo 1. 
Solid line from numerically integrating equation (|8j, dotted lines 
asymptotic limits, equations ((2]l and JSj, short dashed line using 
the original fitting function, equation JgJ, and long dashed line 
the modified fitting function, equation mOI I . Stars from numerical 
simulations. The dot-dashed line is the radius of the mini-halo. 
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Figure 2. The fractional energy input, AE/E, from simulations 
of the best fit Plummer profile for halo 1 as a function of relative 
velocity for encounters with a perturber of mass M* = IMq . The 
open circles have an impact parameter of 6 = 10"'^ pc, while the 
filled circles have b = 10~^ pc. The fractional energy input scales 
as D~^ as expected (lines of gradient —2 have been added to aid 
the eye). 
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In the h ^ R regime the energy input in the simulations 
is significantly larger than that from the analytic impulse 
approximation calculation. This may be indicating that in 
this regime, due to the assymetry of the interaction, the 
(5v).v term in the total energy input does not average to 
zero. It would be interesting to examine whether the energy 
input for b ^ R depends on the mini-halo density profile. 

In Fig. 4 we plot the fractional energy input from an 
interaction with a fiducial star with mass Af* = IMq and 
relative speed v = 300kms~^ for the the best fit profiles 
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Figure 3. The fractional energy input, A.E/ E, from simulations 
of the best fit Plummer profile for halo 1 as a function of perturber 
mass for encounters with a relative velocity v = 100 km s~^ at 
an impact parameter of fe = 10~^ pc. The fractional energy input 
scales as M^ as expected (a line of gradient 2 has been added to 
aid the eye). 
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Figure 4. The fractional energy input in an interaction with 
a fiducial star with mass Mi, = IMq and relative speed v = 
300 km s^^, (Ai?(b)/E)fitj, calculated using equation JSjl for the 
best fit Plummer (solid), CIS (dotted) and NFW (dashed) profiles 
for (from bottom to top) haloes 1, 2 and 3. The dot-dashed lines 
are the radii of the mini-haloes (0.03 pc for halo 1, and 0.008 pc 
for haloes 2/3). 



for all 3 haloes calculated using equation (|8]). The fractional 
energy input for close interactions, which is proportional to 
P^M/ER, varies by a factor of ~ 3 for a given halo and is ~ 
100 times larger for the lighter haloes 2 and 3. This indicates 
that smaller, lighter mini-haloes are far more susceptible to 
disruption by close encounters. For large impact parameter 
interactions (6 S> R) the fractional energy input, which is 
proportional to c?MR^ /E varies only weakly (by a factor of 
~ 3) between haloes and profiles, with the spread in values 
for different profiles for a given halo being comparable to 
that for different haloes for a fixed profile. 

This behaviour can be qualitatively understood by con- 
sidering the asymptotic fractional energy input for a uniform 
density sphere (with p = po for r < R and p = otherwise) : 
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AE{b) 
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1 



{h:s>R), 

(b<^R). 



On galactic scales the red-shift at which a given scale goes 
non-linear, and hence the characteristic density of typical 
haloes, is strongly scale dependent. The comoving scales 
corresponding to the mini-haloes {k > 0(0.1 pc~^)) en- 
tered the horizon during the radiation dominated epoch, 
where (CDM) density perturbations grow only logarithmi- 
cally. The size of the density perturbations, at fixed red-shift, 
on these scales is therefore only logarithmically dependent 
on the comoving wavenumber. Consequently the red-shift 
at which a given physical scale goes non-linear, and hence 
the characteristic density of the resulting haloes, is only 
weakly (roughly logarithmically) dependent on the scale (see 
e.g. Green et al. 2005). Neglecting this weak scale depen- 
dence and making the approximation that po ~ const, then 
AE/E ~ const for 6 > i? and AE/E ~ M"^/^ for b <^ R. 
These scalings are in broad agreement with the trends found 
for the three sample haloes. The weak scale dependence of 
the red-shift of non-linearity will lead to more massive haloes 
typically having lower characteristic densities and hence be- 
ing slightly more susceptible to disruption. This scale de- 
pendence is relatively small however and is comparable in 
magnitude to the dependence on the mini-halo density pro- 
file. 

This behaviour, along with the results from the nu- 
merical simulations for the Plummer sphere, indicates that 
a reasonable approximation to the fractional energy input 
is given by a sudden transition between the asymptotic 
b « / » R regimes: 



V E Ad 

where 
V E /fid.s 



fld.s V b ) 



(AE\ ^Ipc 

V E jfld.o V bs 



W 



G^MlMR^ 



b> bs, 
b<bs, 

1 X 10~ 



(17) 



3 (300kms-^)2(lpc)4 

is the asymptotic large b slope and the transition between 
the two regimes occurs, for the 3 sample haloes, at 

ft, = /-^J R^A^^^R^{0.3-0A5)R. (18) 

The large variation in the value of /3^ for the different 
haloes/profiles has a relatively small (less than a factor of 
2) effect on the value of A, however taking this R/M de- 
pendence into account is crucial for obtaining the correct 
fe <C -R asymptotic behaviour. The NFW profile, however, 
is slightly problematic. As discussed above, its (very large) 
asymptotic value of /3^ is only reached for tiny impact pa- 
rameters. A reasonable prescription for this profile is to use 
the asymptotic fractional energy input at & = O.IR to cal- 
culate the value of /3^ . 

It should be emphasised that the representative sample 
haloes studied in detail by Diemand et al. presumably form 
from 'typical', ~ 1 — 2a, fluctuations. Similarly, our discus- 
sion (above) of the scaling of the energy input with the halo 
mass implicitly assumed that the haloes form from similar 
sized over-densities. Mini-haloes which form from rarer large 
over-densities will be denser and hence more resilient to dis- 



ruption (Berezinsky et al. 2003, 2006; Green et al. 2004, 
2005). More specifically, in the spherical collapse model, a 
halo forming on a given comoving scale from an Na fiuctu- 
ation will have R oc 1/A'', M ~ constant, and characteristic 
density p oc A'^"' (Green et al. 2005). We therefore expect 
that the fractional energy input in close encounters will be 
far smaller for haloes formed from rarer, larger, fluctuations. 
The quantitative effect on the fractional energy input will 
depend on exactly how the characteristic density and den- 
sity profile scale with the size of the overdensity from which 
the mini-halo forms. 



3.4 One-off disruption 

We now use the criterion AE{bc)/E — 1 and the sudden 
transition approximation developed in Section [3.31 above, to 
estimate the critical impact parameter 6c, below which the 
energy input in a single encounter is larger than the binding 
energy. Taken at face value, an energy input AE{b)/E > 1 
might appear to imply that the mini-halo is completely dis- 
rupted. In reality, however, the reaction of a system to a 
sudden change of energy, and in particular the relationship 
between the energy input and the mass lost, is non-trivial 
(see e.g. Aguilar & White 1985; Goodwin 1997; Gnedin & 
Ostriker 1999; and, for the specific case of mini-halo interac- 
tions with stars, Goerdt et al. 2006). The system will expand 
and attempt to revirialise, and during this process two-body 
encounters will redistribute energy between particles. The 
simple criterion AE{bc)/E = 1 allows us to make an esti- 
mate of the impact parameter below which a mini-halo will 
lose a substantial fraction of its mass in a single encounter 
(which for compactness we refer to as 'one-off disruption'). 
A detailed calculation of the mass loss, however, requires 
numerical simulations of the revirilisation and energy redis- 
tribution processes (c.f. Goerdt et al. 2006). 

One-off disruption can not occur if the asymptotic frac- 
tional energy input as b tends to zero is less than one. This 
is the case if 



M* 300 km s" 










Otherwise, 
be _ /A^y/* /M*300kms-^ 

Tpc ^ v~yfld,s Imq V 



1/2 



(19) 



(20) 



In Fig. 5 we plot 6c as a function of M*/w for the best fit 
Plummer profile for halo 1 calculated using the analytic ex- 
pressions derived from the sudden transition approximation 
(eqns. I19ll20p . We also plot the energy input in numerical 
simulations, demonstrating that the sudden transition ap- 
proximation provides a good fit to the transition between 
the {AE{b)/E) > 1 and {AE{b)/E) < 1 regimes. The crit- 
ical impact parameter is quite sensitive to the properties of 
the perturbing star. Thus a full calculation of mini-halo dis- 
ruption will have to take into account the stellar mass and 
velocity distributions. 
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Figure 5. The maximum impact parameter for which one-off 
disruption can occur, be, for the best fit plummer profile for halo 
1 (which has R = 0.03 pc and Vp = 0.013 pc). The solid line 
shows the analytic calculation using the sudden transition ap- 
proximation, equation ||20|I . The symbols are the results of nu- 
merical simulations: open circles for {AE(b)/E) < 0.05, filled 
circles for {AE(b)/E) > 1 (potential one-off disruption) and the 
size of crosses is proportional to {AE(b)/E) in the intermediate 
regime. 



the disc will undergo a change in its energy which is large 
compared to its binding energy and hence lose a significant 
fraction of its mass. The density of stars in the spheroid is 
significantly smaller, ~ W'^Mq pc"'', and declines rapidly 
with increasing Galactocentric radius, therefore mini-haloes 
on orbits which never pass through the disk are extremely 
unlikely to experience a close encounter which removes most 
of their mass. Most mini-haloes will however be on interme- 
diate orbits and spend some fraction of their time passing 
through the disc. For instance a mini-halo on a circular polar 
orbit at the solar radius with speed v — 270 km s^^ would 
spend a fraction ~ 0.08 of its time within the disc, giving a 
disruption timescale of 6 Gyr. Therefore, in the inner regions 
of the MW where orbits pass through the disc, the timescale 
on which 10~^Mq haloes which experience significant mass 
loss in a single interaction is of order the age of the Milky 
Way and a more sophisticated calculation is required. 

Generalising to the more realistic case of a population 
of perturbers with a range of speeds and masses the rate at 
which interactions with impact parameters smaller than the 
critical impact parameter for potential one-off disruption to 
occur becomes 



dJV 
"d7 



dAf*d?; 



■Kvb^{M^/v)dM-^dv, 



(22) 



4 IMPLICATIONS AND OPEN ISSUES 

4.1 Disruption timescales 

As discussed in the introduction (see also Zhao et al. 2005b; 
Goerdt et al. 2006; Angus & Zhao 2006) an accurate calcu- 
lation of the mini-halo survival probability distribution will 
require the combination of simulations of mass loss with or- 
bits in a realistic Galactic potential. In this section though, 
we use the results of our energy input studies in Section 3, 
in particular the sudden transition approximation, to esti- 
mate the disruption timescales for typical mini-haloes as a 
function of mass for some benchmark orbits. 



where d^n/dM^dw is the number density of stars with mass 
between M* and M^-l-dM* and relative speed between v and 
V + dv. We assume that that mass and speed distributions 
are independent so that 



d^n 



dni dn2 



dM^^dv dM^, dv 



(23) 



For the mass distribution we use the Kroupa (2002) stellar 
mass function (MF) 



drii 
dAL 



M, 
M. 



-1.3 
-2.3 



0.08 < M*/Mq < 0.5 , 
0.5 < MJMq < 50, 



4.1.1 One-off disruption 

For the simplified situation where all perturbers have the 
same velocity and mass then the rate at which encounters 
with impact parameter smaller than foe, the critical value for 
which the energy input is larger than the binding energy, 
occur is 



dN 
"dF 



■ Tvnvb^ 



(21) 



where n is the perturber number density. Taking the stel- 
lar mass and relative speed to be fixed at Af* = O.5M0 and 
V = 270kms~^, halo 1 will never undergo one-off disruption, 
while for haloes 2 and 3, using the sudden transition approx- 
imation, the critical impact parameter for one-off disruption 
is 0.0075 pc. Taking a disc mass density of 0.023 Mq pc^'^,[j 
we find a one-off disruption timescale, fdis ~ l/(dA^/di), for 
halo 2/3 of 0.5 Gyr for (rare) halo orbits which lie entirely 
within the Galactic disc for foe — 0.0075 pc. This indicates 
that a 10~^Mq mini-halo which spends most of its time in 



^ This corresponds to a surface density of 46 Mq pc ^ (Kuijken 
& Gilmore 1989) over a height of 2 kpc. 



dn2 
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2cr| 


d« 


- (27r)i/2a. ""P 



which is a good fit to the local field population (see also 
Chabrier 2001). We ignore the contribution from brown 
dwarfs as, due to the M^ factor, this population - whilst nu- 
merous - makes only a small contribution to the disruption 
rate. We normalise the mass function so that the total mass 
density is 0.023 AIq pc"''. We take the speed distribution to 
be Gaussian about the mini-halo speed, V = 270kms~^, 



(24) 



with stellar speed dispersion cr* — 25kms^^. 

The resulting one-off disruption timescales are 0.8 Gyr 
for halo 1 and 0.5 Gyr for halo 2/3. For halo 2/3 the dis- 
ruption time is similar to that calculated assuming delta- 
function mass and speed distributions. The main result 
though is that, once the spread in stellar masses is taken 
into account, the more massive halo 1 can undergo one-off 
disruption on a timescale smaller than the age of the MW. 
Taking into account the spread of masses and velocities is 
therefore crucial for calculating the mass threshold above 
which mini-halos will not lose a significant fraction of their 
mass in a single encounter. 
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4.1.2 Disruption through multiple encounters 

The timescale on which a mini-halo will lose a significant 
fraction of its mass as a result of the cumulative effects of 
encounters with AE{b)/E < 1 (which, for compactness, we 
refer to a 'disruption through multiple encounters') can be 
estimated as 

*di== = .,F/,,^ • (25) 

(d_B/dt)tot 

This is likely to be an overestimate; the mini-halo density 
profile changes in response to interactions and this appears 
to reduce the effect of cumulative interactions (Goerdt et al. 
2006). 

Starting, once again, with the simplifying assumption 
that all stars have the same mass and relative velocity then 



{dE/dt)t 
E 



= 2-K 



nv — -^—^bdb^ 
E 



(26) 



and taking the same parameter values as above we find, for 
orbits which lie entirely within the disc, tdis = 0.4 Gyr for 
halo 1 and idis = 0.5 Gyr for halo 2/3. The shorter timescale 
for multiple disruption for halo 1 reflects the fact that it can 
not undergo one-off disruption and hence 6c = 0, whereas 
for halo 2/3 be = 0.0075 pc. 

Generalising to a distribution of masses and relative 
velocities the fractional energy input rate becomes 



(d£/df)tot 



i)c(M*/t)) 



= 27r 

d^n 
dA/*d« 



vb — T^db 
E 



dudM* 



We now find tdis = 0.6 Gyr for halo 1 and idis 
haloes 2/3. 



(27) 
0.5 , Gyr for 



The net energy input rate is the sum of the energy input 
rates from 'one-off' and 'multiple disruption', and the net 
disruption timescale will be shorter than the characteristic 
timescales for 'one-off' and 'multiple disruption' individu- 
ally. 

We have assumed that the stellar density within the disc 
is uniform. In general, clustering will increase the spread in 
disruption timescales for mini-haloes of a given mass. In ad- 
dition most stars (certainly those with Af* > O.5M0) are in 
fact in binary systems (e.g. Goodwin et al. 2006) and will 
cause a greater disruptive effect than a single star. Systems 
whose separations are significantly less than the mini-halo 
radius (< 1000 AU) will effectively combine the primary 
and secondary masses and, due to the M^ dependence of 
the energy input, even fairly low-mass secondaries may play 
an important role. We estimate that ~ 30 — 40% of stars 
with Mi, > IMq may have a large enough companion to in- 
crease the energy input by a factor > 2|j Even very low mass 
stars (< O.SMq) have a binary frequency of ~ 30% (Fis- 
cher & Marcy 1992), and so the fraction of M-dwarfs with a 
companion that could very significantly increase the energy 



® A companion with a mass ratio greater than 0.4 will increase 
the energy input by more than (1.4)^ ~ 2, and we assume a binary 
fraction of r^ 60% (see Duquennoy & Mayor 1991). 
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Figure 6. The disruption timescales for an orbit entirely within 
the disc as a function of mini-halo mass. Solid line: one-off dis- 
ruption assuming b^ oc M^'^'^, dotted: one-off disruption assum- 
ing 6s oc M^'^ , short dashed: multiple disruption assuming 6s oc 
M^'^^, long dashed: multiple disruption assuming fes ex: M^'-^. 



input is ~ 15 — 20%. Thus, an accurate calculation of mini- 
halo disruption will have to combine simulations of mini-halo 
orbits in a realistic potential with an accurate model of the 
stellar distribution, including the binary fraction, within the 
disk. 



4.2 Mass dependence 

We have seen that more massive mini-haloes are less suscep- 
tible to disruption. It is therefore interesting to investigate 
the mass dependence of the disruption timescales. Further- 
more the WIMP damping scale, and hence the mass of the 
smallest mini-haloes, depends on the properties (elastic scat- 
tering cross-section and mass) of the WIMPs. For generic 
WIMPs Green et al. (2005) found a spread in the minimum 
mass of several orders of magnitude, while Profumo, Sig- 
urdson & Kamionkowski (2006) have recently found that in 
the Minimal Supersymmetric Standard Model the minimum 
mass may vary between 1O~^'^M0 and 1O~*M0. 

The correct mass dependence of the small-b energy in- 
put, and hence the transition radius, 6s, is not known. We 
consider two scalings which should give an indication of the 
trend, and also the uncertainties. Firstly, motivated by the 
the qualitative understanding of the mass dependence of the 
asymptotic limits of the fractional energy input found in sec- 
tion 3.3, we once more consider a uniform density sphere. For 
uniform density spheres (with constant density) the profile 
parameter A{— 4a^/9/3'^) is independent of the radius/mass 
and hence bs (x R oc M ' . For the 3 sample haloes A de- 
creases slightly with increasing mass (albeit with significant 
scatter between profiles and between haloes 2 and 3), and 
6s oc M ■ . We use both of these scalings, normalising in 
both cases to 6s = 0.004 pc at M = 10"'' M©. The correct 
variation might be significantly different from either of these 
scalings, however, and needs to be determined from the pro- 
files of simulated haloes with a range of masses. In Fig. 6 
we plot the resulting disruption times for an orbit entirely 
within the disc (for other orbits the disruption timescales 
scale roughly as the fraction of time spent within the disc) 
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as a function of mini-halo mass, using the sudden transition 
approximation with (A_E/_E)fld,s ~ 10~*. 

For very small mini-haloes, M < 10~'^Mq, the one- 
off and multiple-encounter disruption timescales are in- 
dependent of mass, and roughly equal. The mass inde- 
pendence is because for these small mini-haloes the tran- 
sition impact parameter is smaller than the critical im- 
pact parameter for the range of M* and v values con- 
sidered, bs < bc{M^/v), so that bc{Mi,/v) lies in the 
AE{b)/E oc 6~^ regime and is independent of the mini- 
halo mass. The approximate equality of the one-off and 
multi disruption timescales can be understood by consid- 
ering the simplified case of a delta-function mass/velocity 
distribution once more. Then, using the sudden-transition 
approximation, both disruption timescales are equal to 
[■Knv{M./MQ){300kms-''/v)iAE/Eff/l^{lpcf]-'' if b, < 
6c. The more rapid scaling of 6a, with M (M°-^^ versus M"'^) 
also leads to larger values of idis. The size of these differences 
increases with increasing mini-halo mass. The exact disrup- 
tion timescales of more massive mini-haloes will depend on 
the mass dependence of the impact parameter at which the 
transition between close and distant encounters occurs and 
also how rapidly this transition occurs. 



4.3 Mini-halo radius 

In common with other studies (Zhao et al. 2005ab, Moore 
et al. 2005), we have taken the mini-halo radii to be the ra- 
dius at which the density is 200 times the critical density at 
z = 26 (the red-shift at which Diemand et al. (2005) stopped 
their simulations and plotted the profiles of their sample 
haloes), r2oo(z = 26). The densities of simulated haloes do 
not decline sharply to zero beyond a given radius however 
and, if the mini-haloes remained isolated beyond this red- 
shift their nominal radii (and hence masses and binding en- 
ergies) would increase as the background density decreases. 
As an extreme example, if halo 1 remained isolated to z=0 
then, assuming a NFW density profile, its present day radius 
would be ~ 0.7 pc, its mass would a factor of '--^ 3 larger and, 
using equations ([2| and ((5J, the fractional energy input in 
small (large) b encounters would be substantially decreased 
(increased). The value used for the mini- halo radius there- 
fore has a potentially significant effect on calculations of the 
fractional energy input. 

Once a mini-halo is accreted onto a larger halo it no 
longer accretes further mass onto itself and it is also subject 
to the tidal field of the parent halo. For a mini-halo orbit- 
ing within a Milky Way-like parent halo the tidal radius is 
only comparable to r2oo{z — 26) for very small, of order 
a few kpc, Galactocentric radii. The radius of a mini-halo 
which does not pass through the very central regions of the 
Milky Way will be the smaller of the tidal radius and the 
radius at the time of accretion (both of which are larger 
than r2ooiz = 26)). The red-shift at which accretion oc- 
curs will, however, be different for different mini-halos with 
the same initial properties. A detailed calculation of mini- 
halo evolution will therefore have to include the mini-halo 
merger histories. The majority of mini-halos, in particular 
those which pass close to the solar radius and are hence 
most relevant for WIMP direct and indirect searches, will 



be accreted onto larger halos not long after z=26 and hence 
'"200(2 ~ 26) should be a reasonable estimate of their radii. 



5 DISCUSSION 

We have studied the energy input into earth mass mini- 
haloes in interactions with stars. Using the impulse approx- 
imation (see Spitzer 1958; Gerhard & Fall 1983; Carr & 
Sakellariadou 1999) we have calculated the energy input as 
a function of impact parameter for a range of mini-halo den- 
sity profiles. We also used the DRAGON code (e.g. Goodwin 
et al. 2004a, b; Hubber et al. 2006) to simulate interactions 
with Plummer sphere haloes. We found excellent agreement 
with the impulse approximation in the asymptotic limits 
b <^ / ^ R (where 6 is the impact parameter and -R is 
the mini-halo radius) with a rapid transition at 6 ~ 0.1-R 
between these regimes. We also verified the scaling of the 
fractional energy input with stellar mass and relative veloc- 
ity. 

Using analytic calculations we find that the fractional 
energy input for large impact parameters, b ^ R, appears 
to be fairly independent of the mini-halo mass, varying by 
a factor of -^ 2 for haloes with masses which differ by a 
factor of ~ 50 with a similar variation for different density 
profiles. This behaviour probably reflects the fact that the 
haloes form at roughly the same time and hence have sim- 
ilar characteristic densities. The fractional energy input in 
the 6^0 limit depends quite strongly on the mini-halo 
mass (being larger for lighter haloes) and is also dependent 
on the central density profile. For the NFW profile, which 
has asymptotic inner density profile p oc r~^, the fractional 
energy input only becomes significantly larger than that for 
the cored density profiles at tiny, and hence extremely rare, 
impact parameters i.e. 6 ^ 10"'^ -R. This divergence is there- 
fore essentially unimportant for our calculations, however 
the central regions of haloes with cuspy density profiles may 
be able to survive even after substantial energy input/mass 
loss (e.g. Moore et al. 2005; Goerdt et al. 2006). Motivated 
by the results of our analytic and numerical calculations 
we formulate a fitting function for the fractional energy in- 
put as a function of impact parameter, which we refer to 
as the 'sudden transition' approximation. The slope of the 
fractional energy input at large impact parameters is con- 
stant, while the impact parameter which characterises the 
transition between the limits is mini-halo mass dependent. 

We also investigated the dependence of the critical im- 
pact parameter, 6c, for which the energy input is larger than 
the mini-halo binding energy on the mini-halo mass and also 
the relative speed and mass of the interacting star. As ex- 
pected from the fractional energy input calculations, for slow 
encounters with massive stars 6c is independent of the mini- 
halo mass. There is a critical value of (M*/ii), which in- 
creases with increasing halo mass, below which the energy 
input is always smaller than the binding energy. For all val- 
ues of (M*/M0)(3OOkms~^/«), the results of our Plummer 
sphere simulations are in good agreement with the analytic 
expressions for 6c from the sudden transition approximation. 

We then use the sudden transition approximation to es- 
timate the timescales for one-off and multiple disruption for 
mini-haloes in the MW as a function of mini-halo mass, us- 
ing the approximate destruction criterion AE/E — 1. We 
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take into account the stellar and velocity distribution and 
note that binary stars can cause a significantly greater en- 
ergy input than single stars, due to their greater efi'ective 
mass. For light mini-haloes with M < O(10~'' Mq) the dis- 
ruption timescales are independent of mini-halo mass and, 
for a mini-halo in the inner regions of the MW on a typ- 
ical orbit which spends a few per-cent of its time passing 
through the disc, are comparable to the age of the MW. 
For more massive mini-haloes, M > 0{W~'''Mq), the dis- 
ruption timescale estimates increase rapidly with increasing 
mass, suggesting that the majority of these mini-haloes will 
not be disrupted by stellar encounters. It is important to 
caution, however, that the relationship between the energy 
input and the change in the bound mass is not straight for- 
ward. In particular the mini-halo density profile evolves so 
that successive multiple encounters are less efi'ective than 
would naively be expected and even if the energy input in 
a single encounter is much larger than the binding energy 
a small fraction of the mass can remain bound (Goerdt et 
al. 2006). Therefore these simple estimates are likely to be 
overestimates of the actual disruption timescales. 

Finally we discussed the dependence of the fractional 
energy input on the mini-halo radius assumed. To be con- 
sistent with other studies (Zhao et al. 2005ab, Moore et al. 
2005), we took the radius to be the radius at which the den- 
sity is 200 times the critical density at 2 = 26, the red-shift 
at which Diemand et al. (2005) stopped their simulations 
and plotted the profiles of their sample haloes. This is a 
somewhat arbitrary definition however; the densities of sim- 
ulated haloes do not decline to zero beyond this radius and 
as the background density decreases the nominal radius in- 
creases. The physical extent/radius will in fact be that at 
the time of accretion onto a larger halo, or the tidal radius 
if this is smaller. The tidal radius within a Milky Way-like 
parent halo is only smaller than r2oo{z = 26), at small Galac- 
tocentric radii, however the majority of mini-halos will be 
accreted onto larger halos shortly after this red-shift, so in 
practice r2oo{z = 26) should be a reasonable estimate of the 
radius of most mini-halos. 

A complete calculation of the disruption of mini-halos 
will need to take into account their merger histories, simul- 
taneously and consistently incorporate disruption due to en- 
counters with stars and tidal stripping. Mini-halos formed 
from rare, large density fiuctuations, will be denser, and 
hence more resilient to disruption, than typical mini-halos 
and this will also need to be included. 
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